Resonant light delay in GaN with ballistic and diffusive propagation 
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We report on a strong delay in light propagation through bulk GaN, detected by time-of- flight 
spectroscopy The delay increases resonantly as the photon energy approaches the energy of a 
neutral-donor bound exciton (BX), resulting in a velocity of light as low as 2100 km/s. In the close 
vicinity of the BX resonance, the transmitted light contains both ballistic and diffusive components. 
This phenomenon is quantitatively explained in terms of optical dispersion in a medium where 
resonant light scattering by the BX resonance takes place in addition to the polariton propagation. 

PACS numbers: 78.20.e,78.47.p 
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Production of 'slow light' is a subject of extensive stud- 
ies nowadays [1] . Propagation of an electromagnetic wave 
through a dielectric medium can be considerably slowed 
down via ballistic or diffusive propagation, due to either 
peculiarities of dispersion or multiple scattering acts, re- 
spectively. The ballistic light propagation through the 
dispersive medium with an absorption resonance exhibits 
strong retardation of the energy transport in the vicinity 
of the resonance [2|, so that the group velocity Vg can 
be formally zero or even negative [S]. In semiconductors 
with exciton resonances, this effect is closely related to 
the exciton-polariton formation |4|. The pulse propaga- 
tion time agrees with that calculated from the polari- 
ton group velocity in different semiconductor materials, 
namely, in GaP:N around an isolated bound exciton (BX) 
line [5[ , CuCl near a free exciton resonance ^] , and ZnO 
in the BX vicinity . 

Light diffusion involving multiple photon scattering 
has been observed in various disordered media, including 
dielectric micro-spheres, polymers, and ZnS nanocrystals 
[8|. It is accompanied by intriguing phenomena, such as 
weak localization and coherent backscattering of light [9| . 
A strong slow-down of light transfer due to diffusive mo- 
tion of photons was discussed in, e.g., [10] . Resonant 
back-scattering of light has been observed in quantum 
well structures [ll| . On the other hand, in bulk semicon- 
ductors, no experimental evidence for slow-light forma- 
tion due to multiple scattering has been reported so far, 
to our knowledge. 

In this Letter we report the first observation of a strong 
delay of light propagation in a high quality bulk GaN 
crystal, which amounts to 470 ps for a 1 mm long way, 
corresponding to a light velocity of about 2100 km/s. We 
demonstrate a resonant nature of this phenomenon and 
show that both light diffusion and retarded ballistic prop- 
agation of exciton-polaritons can contribute to the delay. 
Due to the wide use of GaN in modern optoelectronics. 
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FIG. 1: (Color onhne) (a) Backward and (b) forward TR 
PL images recorded in the 1-mm sample with excitation by 
a 4.66 eV laser line. In the colored images, the PL intensity 
(logarithm scale) increases with the color variation from dark- 
blue to red. The inset presents time-integrated spectra of 
backward PL in the 1-mm sample (1) and forward PL in the 
1-mm (2) and 2-mm (3) samples. 



this effect is promising for practical applications. 

The studies have been done by time-of-flight spec- 
troscopy using free-standing GaN epilayers grown by 
HVPE, which are among the best available. They exhibit 
photoluminescence (PL) spectra with the width of neu- 
tral donor BX lines about 0.5 meV. We focused on three 
samples of different thicknesses, L\ 1, 2 and 0.3 mm. 
The first sample has the carrier density of 8x10-*^^ cm~^, 
two others around 5x10^^ cm~^. In these experiments, a 
Hamamatsu streak camera with a ~20 ps temporal res- 
olution was exploited to record the time-resolved (TR) 
images. The third and second harmonics from a tuned 
Ti: sapphire femtosecond pulsed laser were used as a light 
source. The temporal width of the detected laser pulse 
(20-50 ps) was determined by the instrumental accuracy. 
The transmitted light was registered in the forward di- 
rection along the normal to the surface. Conventional 
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FIG. 2: (Color online) TR images (logarithm scale) of 3.486 
eV pulses recorded (a) without a sample, (b) through the 1- 
mm sample, (c) The pulse shape at 3.465 eV calculated for 
ballistic (1) and diffusive propagations in the sample neglect- 
ing (2) and taking into account (3) the 97% light reflection 
at boundaries, (d)-(f) TR images of 3.464 eV pulses passing 
through the sample at different temperatures. 
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FIG. 3: Normalized angular diagrams of the light intensity 
measured at different photon energies in the 1-mm sample 
(and without that). The inset presents deconvolution of a 
diagram into narrow and broad components. The width of 
the ballistic part is determined by the light scattering at the 
sample surface. Some distortion of the broad component at 
large angles is due to the shading by cryostat windows. 



backward PL was measured for reference. 

With excitation in a highly absorbing region by a 4.66 
eV laser line, the backward TR PL images (Fig. 1 a) 
show a fine excitonic structure with the pronounced A 
and B free exciton lines and the BX doublet at 3.472 
eV, which is due to the exciton localization at Si and O 
donors. The transmitted PL images turn to be absolutely 
different. They exhibit a sharp quenching of the signal 
at a boundary located just below this doublet (Fig. 1 b). 
This boundary at ~3.47 eV is rather sharp in the best 
1-mm sample; in other samples it is more diffuse and red- 
shifted. The other remarkable feature is the bending of 
the PL streak towards this boundary. It means that the 
emission leaving the sample appears with a delay in time, 
which increases with approaching the BX resonances. 

When the laser pulse was tuned to 3.486 eV to be res- 
onant with the free excitons, the transmitted signal in- 
cludes the excited PL along with the non-absorbed part 
of the pulse. Both components are cut and bent in a sim- 
ilar way (Fig. 2, b). The excited PL was negligible, when 
the excitation energy was tuned to 3.464 eV, i. e. be- 
low the exciton resonances. This permits us to study the 
propagation of the laser pulse solely. Again, we observe 
the similar bending towards BX. The images recorded at 
different temperatures show a shift of the cut-off with 
the increase of temperature. This shift is ~10 meV in 
the 2 — 70 K range, which is consistent with the variation 
of the band gap in GaN and, hence, of the BX energy. 

There is one more important feature - an additional 
curved streak ("reflex"), which is clearly seen in the im- 
ages without PL taken with high resolution (Fig. 2 d-f). 
Its bending corresponds to the threefold delay time of the 
first transmitted pulse. We assume that this feature ap- 
pears due to the triple coherent passage of a wave packet 
(forward-back-forward) provided by the light reflection at 
the sample boundaries. This seems to indicate that the 



pulse propagation in our samples is predominantly ballis- 
tic, otherwise a photon should lose the direction memory 
after a few acts of scattering. The diffusion model pre- 
dicts an extension of the pulse in the time domain but 
not an extra peak. This is illustrated by Fig. 2 (c), which 
presents the calculated pulse shapes, based on the models 
described below. 

On the other hand, the reflex is clearly seen only in 
the best 1-mm sample. In other samples with higher im- 
purity concentration, it is less pronounced; usually only 
a smooth cloud is recorded. Besides, the reflex is not 
visible close to BX and disappears with the temperature 
rise, while the delay still exists. All this shows that a 
non-ballistic contribution to the delay of light in the BX 
vicinity cannot be excluded. 

In an attempt to gain further information, we have 
measured the far-field angular dependence of the trans- 
mitted light. We expect the ballistic propagation to yield 
a narrow angular diagram of the output intensity, while 
the diffusion is likely to provide a much broader distri- 
bution. For this experiment, we used a Xe lamp, whose 
light passed through a monochromator impinging nor- 
mally to the sample surface. The output intensity was 
registered at different angles by means of rotation of a 
stage, holding the lamp and a cryostat, with respect to 
a detecting system. As one can see from Fig. 3, the an- 
gular distribution of the light intensity clearly indicates 
a superposition of two components - relatively narrow 
and broad ones. We ascribe them to the ballistic and 
diffusive parts of the transmitted light, respectively. We 
found that the intensity ratio between two components 
depends on the energy. The ballistic part dominates in 
the region of relative transparency, while the diffusive 
mechanism plays the major role in the vicinity of the BX 
resonances. 

The coexistence of two regimes of light propagation is 
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FIG. 4: Delay times measured in the 1-mm (squares), 2- 
mm (stars), and 0.3-mm (triangles) samples. Open squares 
are the data on the reflex in the 1-mm sample. The fitting 
curves are presented for ballistic (solid lines) and diffusive 
(dotted lines) propagation. The inset presents reflection (R) 
and transmission (T) spectra, measured using a lamp in the 
1-mm sample, with their simulations (dotted lines). 



cooj equal to 3.4720, 3.4785, 3.4837, and 3.5016 eV for 
BX, A, B, and C excitons, respectively. For either of the 
resonance, Tj is taken as 0.012 meV, which is of the order 
of the homogeneous line width of an exciton-polariton in 
other semiconductors [l3[ . The inhomogeneous width Aj 
is taken as 0.75 meV, close to the average PL linewidth 
in our spectra. For the A, B and C excitons, / equals to 
0.0026, 0.0017, and 0.00037, respectively, which is com- 
parable with the previously reported values for GaN [l2| . 
The oscillator strength for the bound exciton with the 
donor concentration N can be estimated as [l4 
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where mo is the free electron mass, e is the elementary 
charge, Pcv is the interband matrix element, and "^{r, p) 
is the envelope of the electron-hole pair function depend- 
ing on the center of mass position r and the relative 
electron-hole vector p. Under the assumption that the 
exciton is localized as a whole it has the form 



a surprising observation. It is well-known that the bal- 
listic propagation time Tb{oo) depends linearly on the 
sample thickness Tb{uj) = L/vg{uj)^ while the delay 
time caused by a conventional diffusion process would 
increase quadratically with the sample thickness. Thus 
the presence of the diffusive component could result in 
a marked extension of the streaks in the time-resolved 
images, which is not observed. The question arises why 
so different processes yield similar propagation times. To 
clarify that we model both processes separately. 

For the ballistic propagation, the group velocity 
Vg{u;) = duo/dk^ where the wave vector kiuo) = 
{uo/c)^Je{uo). The frequency dependent dielectric con- 
stant in the vicinity of exciton resonances can be written 
neglecting spatial dispersion as 
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Here, £5 is the background dielectric constant, j = 
A^B^C^ BX denotes an exciton resonance. For simplicity 
we consider only one bound state, making no difference 
between an exciton localized at O or Si donors. Each 
resonance is described by a frequency cc;o,j, an oscillator 
strength fj , and a damping term Vj . The inhomogeneous 
broadening of the exciton line is taken into account by 
means of convolution with the Gaussian with the width 
Aj. We underline that only this formalism has allowed 
us to fit the reflection and transmission spectra simul- 
taneously (see inset in Fig. 4). Note that the conven- 
tional approach [l2[ yields a reasonable fitting of reflec- 
tion spectra but strongly overestimates the absorption of 
light near the exciton resonances. 

We have performed a modeling assuming £5 = 9.5 and 
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Here we took the exciton localization radius R = as^ 
where as = 2.8 nm is the Bohr radius, which allows us 
to reproduce the observed value of the BX binding energy 
|l5| . N is taken equal to the carrier density; pcv = 10~^^ 
g cm/s. The best fit is obtained with f ^ 2 x 10~^ 
for the 1-mm sample, what is somewhat smaller than 
an estimation according to Eq. ([2|). For other samples 
f is taken an order of magnitude higher following the 
difference in the carrier densities. Figure 4 demonstrates 
a reasonable fit with these parameters and clarifies the 
origin of the sharp cut-off. The calculated dependencies 
have a strong dip near the BX resonance. It appears 
because the phase and group velocity of light strongly 
decrease here, while the absorption is enhanced [3]. 

To consider the diffusive propagation we assume that 
elastic neutral donor scattering dominates the delay, ne- 
glecting any inelastic scattering by impurities or acoustic 
phonons [Ifi]. We write the diffusion equation for the 
concentration of photons n(cj, x, t) at frequency cj, coor- 
dinate X, and time t 
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Here frequency dependent D{co) and r(a;) are a diffusion 
coefficient and photon life- time, respectively. The life- 
time is defined as r{uj) = tq • (/^(cc;)//o), where Iq and tq 
are a photon mean free path and time between two acts of 
scattering in the region of relative transparency; Ia{^) is 
the absorption length. All the photons are assumed to be 
initially created in a thin layer near the sample surface: 
n(c(;,x,0 = n{u;^0)S{x). We have performed calculations 



of the pulse propagation in the diffusive model neglecting 
and taking into account the light reflection at boundaries 
and found no qualitative difference (Fig. 2, c). Below 
the light reflection is neglected. The delay time due to 
diffusive photon motion, To (a;), can be found from the 
condition dn{u; ^ L ^ t) / dt = which corresponds to the 
intensity maximum. In the general case it is written as 
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We have the two following limiting cases: i) small absorp- 
tion, t{uj) ^ L'^/D{uj) where the delay scales as a second 
power of the sample thickness, Td(cj) = L'^ /[2D{u;)]; ii) 
strong absorption, r{u;) <C L'^ /D{uo)^ where the delay is 
proportional to the sample thickness as with the ballis- 
tic propagation, Td{uj) = {L /2) ^ t{(jo) / D{(jo) . The latter 
result is a consequence of the competition between the 
'fast' absorption and 'slow' diffusion processes. 

The diffusion constant D{(jo) = Vg{uj)ltr{^)/^, where 
ltr{^) = {^^{^))~^ is the photon mean free path deter- 
mined by the light scattering cross-section cr(a;) and the 
concentration of scattering centers N. Assuming that the 
light scattering is caused by the bound exciton resonances 
one may write a{u) in a general form [l7[ [l8[ 
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Here, 7 is an effective broadening parameter. The main 
contribution to the delay time comes from the strong fre- 
quency dependence of the cross-section. Strictly speak- 
ing, Eq. (|6]) describes scattering caused by individual 
BX. This is valid only for the diluted donor concentra- 
tion satisfying the condition N'^^^X <C 1, where A is a 
wavelength inside the crystal. Near BX, this is fulfilled 
for N ~10^^ cm~^, while the donor density is higher in 
our samples. On the other hand, the spatial distribution 
of donors is inevitably inhomogeneous and the clusters 
of closely lying BX act as virtual scattering centers with 
effective concentration and localization radius. For the 
sake of demonstration, we present here Toioo) calculated 
as for the individual scatterers. We use 7 = 0.75 meV, 
T^/h = 2xl0~^ s/cm, and / calculated according to Eqs. 
(O, ©. While for the 1-mm sample some agreement is 
achievable with N = S x 10^^ cm~^, for other samples 
we are impelled to take the effective A^=2xl0^^ cm~^ to 
obtain a reasonable fitting (Fig. 4). 

Our simulations of the experimental data are done 
in the limiting cases of either ballistic or diffusive light 
propagation. We abstain from the fitting in a combined 
model; it hardly can be representative, because the ratio 
between these components is too dependent on the sam- 
ple quality. Our main goal is to demonstrate that the 
diffusive mechanism can give the delay of the same order 
of magnitude as the ballistic one, and that its dependence 
of the delay on the thickness can be linear as well. 



In conclusion, we reported here a strong delay in light 
transfer through high quality bulk GaN samples and 
present a quantitative description of this phenomenon in 
the frameworks of optical dispersion and resonant elastic 
scattering models. We believe that our findings have a 
particular importance for the development of a new gen- 
eration of smart optoelectronic devices, which rely on a 
controlled delay of light propagation. 
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